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Abstract 

The radiative pion capture process in nuclei is approached by using a continuum 
shell-model description of the nucleus, together with a phenomenological treatment 
of the two particle-two hole effects. It is found that these effects play an impor- 
tant role to reproduce the observed experimental photon energy distribution. This 
distribution as well as the integrated one depends significantly on the details of the 
mean field potential. This makes this process interesting to investigate the nuclear 
structure dynamics. 
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1 Introduction 



The radiative pion capture (RPC) in pionic atoms has been extensively studied in the 
past In this process, the vr" sited in a given nl orbit of a pionic atom is captured 

by the nucleus emitting a photon. The interest of this reaction is threefold. First, it 
provides direct information about the vr-nucleus interaction, responsible for the formation 
of the pionic bound state. Second, it proves the nuclear structure, specifically the details 
of the proton-distribution, and the nature of the two particle-two hole {2p2h) effects in 
the nuclear medium. Finally, it gives valuable information on the tt^p '-/n reaction on 
the nucleon at very low pion kinetic energies and in this sense it is complementary to the 
pion photoproduction reactions. In this paper we will concentrate on the second aspect 
mentioned above. 

Technically, the most difficult aspect in the theoretical side is the evaluation of the 
sum over final nuclear states. This has made necessary the use of different approaches 
such us the closure approximation [^] or the RPA sum rules [Q. In the first case, the 
results show a considerable dependence on the average photon energy used, a problem 
which is minimized by the RPA sum rules method. However, this two approximations 
only give integrated quantities (e.g. the radiative capture widths and branching ratios), 
but observables such as the photon energy distribution are out of their possibilities. 

More recently, Chiang et al. have proposed a many-body approach in which both 
the photon energy distribution and the radiative capture width can be calculated. Starting 
from techniques developed for the (7, tt) reaction 0, this method permits the use of 
accurate pionic wave functions and the incorporation of the medium polarization effects. 
Besides, the calculations are considerably simplified because they are done in infinite 
nuclear matter and applied to finite nuclei by means of the local density approximation 
(LDA). Though the obtained branching ratios and widths are in overall agreement with 
the experimental results, the photon energy distribution does not show the marked tail 
of the data at low photon energies. This last problem is partially due to the limited 
momentum components available to the nuclear wave functions obtained from the Fermi 
gas (FG) model. 

Shell-model (SM) calculations do not present this deficiency. This was observed in 
the pioneering work of Vergados |l| and Eramzhyam et al. and in the results obtained 
by Ohtsuka and Ohtsubo and Wiinsch |]IU| within the continuum SM approach. Much 
of this work has been devoted to the analysis of the role played by the resonances in the 
reaction and the main conclusion is that the RPC can be described assuming a dominant 
resonance character for the process. 

In this work we want to investigate the dynamics of the RPC reaction in more detail 
by analyzing it with a continuum SM recently developed to study the quasi-free inclusive 
electron scattering by nuclei [|Il|, |I2|. In particular we want to address different questions 



related to the importance of the parameters of the SM potential, to the problem of the 
orthogonality of the single-particle wave functions and to the effect of the inclusion of 
the particle-hole [ph) self-energy in the medium on the ph propagator. This latter effect 
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accounts for the final state interactions (FSI) of the emitted neutron with the residual 
nucleus. 

The organization of the paper is as follows. In sect. 2 we give the details of the model 
used to perform the calculations with the explicit expressions of the necessary matrix 
elements placed in the appendix. In sect. 3 we show the results and comment the main 
features. In sect. 4 we summarize our work and draw our conclusions. 



2 Details of the model 

2.1 The RPC rate 

In this work we study the following reaction: 

where the pion is captured by the nucleus emitting a neutron and a photon, being the 
latter one the only particle detected. 

The unpolarized differential photon distribution for the capture of a pion, bound in 
the orbit nl, from a spin zero nucleus is given by p!3| : 



1 

fJ 



(q) = ^j^i:S{q + Ef-E,)[\{f\\Tf{qm\' + | (/||Tr^'(?) |K) . (2) 



Here it has been taken into account the fact that the nuclear states as well as that of the 
pion have good angular momenta, what makes useful to perform a multipole expansion 
of the operators involved. In Eq. (^, a is the fine structure constant, q is the energy of 
the emitted photon, \i) and |/) represent the initial (taking into account the pion degrees 
of freedom) and final hadronic states of the system and Eij are their respective energies. 
Besides, Tj^ and Tj"^ are the electric and magnetic multipole operators which are given 
by: 

TjUq) = ^ /dV {V X [jj{qr)Yfj{r)]} ■ J(r), (3) 
q J 

TT^Fiq) = /dVj,(gr)Yf,(r)-J(r), (4) 

where jj is a spherical Bessel function and Yj^ a vector spherical harmonic. Finally, J(r) 
is the electromagnetic current operator. 

The amplitude for the elementary process n^p 772, which determines the hadronic 
current, depends [|I^] on the spin of the nucleon involved in the process, the photon 
polarization and the pion and photon momenta. However, the average momentum of the 
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bound pion in light and medium nuclei is small and in this situation it is plausible to take 
into account only the so-called KroU-Ruderman term (which is the one surviving in the 
limit of zero pion momentum). Then we have for the electromagnetic current operator 

J(r) = zv^-£ (l + ^) E^-(O0-(r)^W5(r-r,). (5) 

Here, the index i runs over the A nucleons of the system with coordinates r,, / is the 
pion-nucleon coupling constant (/^/47r = 0.08), is the pion mass, M is the nucleon 
mass and cr are spin Pauli matrices. Besides, the isospin operator t_ transforms a proton 
into a neutron. Finally, 0-(r) is the pion field operator which can be developed in a 
complete basis as follows: 

= E / Mr) Ck, (6) 

k pEl 

where k runs over all the pionic atom orbits, 0^ is the corresponding pion wave function, 
Ck is the pion annihilation operator, and E^. is the pion energy including its mass. 

After substituting the current (|]) in Eqs. @) and (||) it is possible to obtain the tensor 
form of these electromagnetic operators and calculate the matrix elements involved in 
Eq. (0). These expressions can be found in the Appendix. 



2.2 The model of nuclear structure 

The matrix elements in Eq. (j^) require the definition of the initial and final states of 
the system. We are dealing with closed-shell nuclei and then our initial state is built by 
coupling the 7r~ atomic wave function with the closed-shell core. Thus the initial state is 
given by: 

K) = |J,M,) = [|7r-;n/)® 10+)];'' (7) 

with Ji = I and |0^) the Slater determinant corresponding to the closed-shell nucleus. The 
Slater determinant is built with single-particle wave functions obtained by diagonalizing 
a Woods-Saxon potential well 

Vwoods-saxon(r) = K)/(r, Rq, oq) + Vls 1 h Vc{r), (8) 

r ar 

where 

fir,R,a)= \ (9) 

and Vc{r) is the Coulomb potential created by an homogeneous charge distribution of 
radius Rc- The parameters of this potential are adjusted to reproduce different nuclear 
properties such as, e.g., the experimental single-particle energies around the Fermi level, 
the r.m.s. radius or the charge density |T^. On the other hand the pion state, |7r~; nlrrii), 
is defined as follows: 

0_(r)|7r-;n/m;) = -^Rl^{r)Yim,ir)m. (10) 
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with |0)7r the pion vacuum state, n,l and m the principal, orbital angular momentum 
and magnetic quantum numbers of the pion atomic orbit and Yim{r) spherical harmonic 
functions. The radial pion wave function Rnii'f') and the pion energy E^i are obtained 
by solving the Klein-Gordon equation with a pion-nucleus optical potential. This poten- 
tial has been developed microscopically and it is exposed in detail in Ref. for pionic 
atoms. It contains the coulombian 7t~- nucleus interaction and the ordinary lowest or- 
der optical potential constructed from the s- and p-wave vrA^ amplitudes. In addition 
second order terms in both s- and p-waves, responsible for pion absorption, are also 
considered. Standard corrections to both the electromagnetic (finite size of the nucleus 
and vacuum polarization) and strong (second-order Pauli re-scattering term, ATT term, 
Lorentz-Lorenz effect, long and short range nuclear correlations) parts of the potential are 
also taken into account. This theoretical potential reproduces fairly well the data of pio- 
nic atoms (binding energies and strong absorption widths) [|16| and low energy 7r-nucleus 
scattering [0 . 

In Eq. (0), the initial energy Ei corresponds to the energy of the (7rX)bound system, 
including the pion binding energy. As we are only interested in the excitation nuclear 
energy, we take the energy of the initial closed-shell nucleus as our zero of energies. Thus 
the initial energy Ei can be identified as the pion mass plus its binding energy, E^i. 

The final states are described in terms of a neutron particle in the continuum and a 
proton hole coupled to the final angular momentum Jf, 

\f)^\ph-';JjMf)=[al®bi]l'^^ |0+) (11) 

where aj, {b\) is the particle (hole) creation operator. The particle states are obtained by 
solving the Schrodinger equation for positive energies. We have used the same Woods- 
Saxon potential as for the discrete single-particle levels in order to ensure the orthogo- 
nality of the states involved in the process (see Ref. [jlll for details). It is important to 



mention here that, in view of the definition of the final states we are using, the sum over 
final states |/) in Eq. (||) should be changed to: 

E — EE- (12) 

/ Ph Jf 



Therefore, to evaluate Eq. (Q) we must perform a sum over the quantum numbers of 
the final particle-hole states {h = rih, lh,jh, P = ^p, Ipijp and the total angular momentum 
Jf of the ph excitation) and the photon angular momentum (J). The sum over ep is 
indeed an integration which is performed by means of the delta of energies which appears 
in Eq. (|^). Actually, in our model, the final energy Ef is given by 

Ef = Am + ep - en, (13) 

with ea the energy of the single-level a and Am the difference between the masses of the 
particle (neutron) and hole (proton) nucleons. Thus, for each value of the photon energy 
q and each hole h, ep is univocally determined by the delta function in Eq. (^. 
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Nucleus 


Shell 


eh{th) 


(exp) 




lP3/2 


-18.05 


-15.96 


160 


lPl/2 


-12.74 


-12.12 




lP3/2 


-16.92 


-18.45 




2Si/2 


-10.30 


-10.85 




14/2 


-8.71 


-8.33 




14/2 


-16.19 


-12.27 



Table 1: Proton hole single-particle energies (MeV). We show the theoretical energies com- 
puted with the SM potential described in sect. |2]^ (third column) and the experimental 
ones (fourth column). For more details, see Ref. [p!5[ . 

Note also that the sum over the photon multipoles, J, ranges from one to infinity. 
The larger is the photon energy or the nuclear mass, the bigger is the number of photon 
multipoles needed to obtain convergence. In the RPC process, the maximum photon 
energy is around the pion mass and we have found in ^°Ca a good convergence with four 
or five multipoles. 

One final detail concerns to the hole-energies used in our calculation. As it was 
mentioned above, our mean-field potential has been adjusted to reproduce single-particle 
energies around the Fermi level. To improve the energetic balance in our model, in the 
actual calculation the experimental proton single-particle energies values, when known, 
have been used. In Table |l|, we compare these values to the computed mean-field proton 
hole energies. 



2.3 Effects of the 2p2h correlations. 



In sect. |2.2|, we have outlined our continuum SM approach to the nuclear structure. 



In this framework, the final state interactions (FSI) of the outgoing particle (neutron) 
with the residual nucleus are partially taken into account. However, 2p2h correlations 
are a fundamental ingredient to achieve a reasonable descriptions of these FSI. Here, we 
discuss briefly the phenomenological model for such effects used in this work. We follow 
the approach of Ref. |18[. There, the propagation of the mean field ph excitations in 
the nuclear medium is modified by including an approximated complex ph self-energy. 
This self-energy is determined by the effective nucleon-nucleon residual interaction in the 
medium not included in the mean field approach. The inclusion of the 2p2h correlations 
improves our treatment of the FSI and it also accounts for corrections to the pure one- 
body RPC due to two nucleons mechanisms. This approach has proved to be successful in 



the study of the transverse |12| and longitudinal [|I9| response functions in the quasielastic 
region for inclusive (e, e') processes. Also a recent approach |]20|, where the propagation of 
ph excitations in the medium is described in terms of particle and hole spectral functions 
gives similar results to those of Refs. |T^ and |19[|. 
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In previous approaches |l| , FSI were taken into account by putting the outgoing neu- 
tron in a wave of an neutron-nucleus optical potential. The imaginary part of the optical 
potential reduces the outgoing neutron flux and therefore reduces the integrated RFC 
rate. However, we believe this procedure is incomplete: the imaginary part of the opti- 
cal potential is mainly due to two-body mechanisms which also contribute to the RFC 
rate and have to be considered in order to be consistent (see Ref. for more details). 
Our approach, where we incorporate a ph self-energy, automatically takes into account 
both physical processes mentioned above, leading to a redistribution of the strength of 
the differential photon rate but as we shall see leaving the integrated RFC rate almost 
unchanged. 



Following the steps of Ref. to include the 2p2h correlation effects to the RFC 
differential rate, one should fold the pure one body SM differential rate, (^dT^jJi\q) / dq^ , 
given in Eq. (0), with the electromagnetic current and final states defined in Eqs.(BD and 



[TlD , with the imaginary part of the ph propagator: 







i) ^' / [p{m^-q\m^-q)+ p{m^-q',q-m^)] {U) 



dq Jo \ dq' 



with 



1 T^i^uj^ 

^^^'""^ = 2;^(E-c.-AH)2 + (riH/2)2- ^^^^ 

The functions and F''^(cc;) define the complex self-energy, 

EM = AM-^FHa.) (16) 



of a ph with excitation energy lu. The self-energy Ti(uj) has been taken from Ref. |19 



3 Results 

In this section we apply our model to calculate the RFC differential rates for three closed- 
shell medium nuclei, ^^C, ^^O and ^°Ca, and compare our predictions with the available 
experimental data. 

Experimentally, it is rather difficult to distinguish between RFC processes from dif- 
ferent pionic atom orbits. Indeed, only the weighted ratio 

dB^^^u^dT^ 
dq ^T-f dq ^ ' 

can be measured. Here, Uni gives the absorption probability from each nl pionic level, 
taken into account the electromagnetic transitions and the strong absorption, and are 
normalized to the unity. F^^'^ is the total pion absorption width from the orbit nl and 
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Nucleus 


nl 




r»f (KeV) 




Is 


0.1 


3.14 




2p 


0.9 


0.00136 


160 


Is 


0.1 


7.92 




2p 


0.9 


0.00676 




2p 


0.7 


1.59 




3d 


0.3 


0.0007 



Table 2: Values of ujni § and T^^" PI— ^ used in this work. 



V^^i is the width due to the radiative capture of the pion from the orbit nl. As disccused 
in [0, we make the approximation of putting all the weight of radiative capture in the two 
pionic levels observed with the X-ray technique. The values used in this work for these 
parameters are shown in Table |[ In what follows we will refer to dR^'^^ / dq instead of the 
single differential decay rate of Eq. (|^). 

We start now discussing our results. Firstly, we focus on the purely one-body process 
and later we will introduce the effects of the 2p2h correlations. In Fig. |I|, we show the 
contribution to dR^'^''> /dq from the different holes configurations in ^^C, ^^O and ^°Ca. 
The most important contributions are due to the Ip shell in ^^C and ^^O and to the 
Id and 2s shells in ^°Ca. In all cases they correspond to the outest hole configurations, 
which maximize the overlap with the pion wave function. The end of the contribution of 
a particular shell produces a protuberance in the total sum. On the other hand, in the 
plots also appear resonance structures, specially in ^°Ca where turn out to be extremely 
narrow with typical widths of the order of 10 KeV. This was already pointed out in [^. 
However, the resonance structure for low nuclear excitation energies depends strongly on 
the details of the potential. The use of a different nuclear structure model will produce 
drastic changes of both their widths and positions. 

Note also that in the data, these resonances will be convoluted with the experimental 
photon energy resolution, which is much greater than their widths and therefore the sharp 
resonances shapes will smear out and a smoother spectrum will be observed. Here, we are 
mainly interested in the global shape of the photon distribution and the study of these 
fine details is out of the scope of this work, because it will require both a more refined 
nuclear model and more precise experimental data. 

In Fig. 1^ we compare the total sum of all hole contributions shown in Fig. |l| with the 
experimental data. In calcium and oxigen, following the discussion above, we have elimi- 
nated the sharp resonance structures in our differential rates by convoluting the responses 
with a Gaussian weight function /(cu) oc ea;p(— cj^/F^) of width F = l.SMeV. This convo- 



lution procedure keeps unchanged the integrated RFC rate ||2^. The experimental data 
are given in arbitrary units, thus to compare with the experiment we have normalized our 
results to the data at g = 105 MeV. The reason to normalize our results in this region 
of intermediate energies is twofold. Firstly because is far from the low nuclear excitation 
energy region (high photon energies) where, as we discuss above, the theoretical results 
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q [MeV] 

Figure 1: Contribution to dR^'^'> /dq from the different holes configurations in ^^C, ^^O 
and ^°Ca. In all figures, the sohd lines stand for the total sum of the contribution of 
the hole configurations. Each hole can contribute up to a maximum photon energy, 
Q-max = Ei — Am + eh, which is different for each shell. This fact can be used to identify 
the contribution of each hole in the figures: the deeper is the shell, the smaller is ql^ax- 
Prom the deepest to the outest shell, the hole configurations involved are: lsi/2, lp3/2, 
lpi/2, 14/2, 2si/2 and 14/2. 
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Potential Vq (protons) Vq (neutrons) R^'^^[% 



A 


-57.5 


-55.0 


1.25 


B 


-50.0 


-55.0 


1.54 


C 


-50.0 


-50.0 


1.87 


D 


-50.0 


-60.0 


1.21 



Table 3: Different Woods-Saxon potentials considered in this work. The second and 
third columns give proton and neutron Wood-Saxon well depth parameters (in MeV) for 
^"Ca used to compute the different curves appearing in Fig. ^. The rest of the parameters, 
not shown here, are common for all potentials and can be found in Ref. |T^. The po- 
tential B has been adjusted to reproduce the experimental proton density, whereas C(D) 
corresponds to a 10% decrease (increase) of the depth for the neutron Wood-Saxon well. 
The fourth column gives the integrated ratio R^'^^ (in percent) for each potential. 

depend strongly on the nuclear model. Secondly, because the low photon energy region is 
sensitive to the high momentum components of the hole wave function, and one expects 
theoretical uncertainties in the nuclear model. Furthermore, as we will show below, the 
effect of 2p2h correlations should be also included in the low energy photon region. 

As it can be seen in Fig. ^ our description gives the gross features of the experimental 
data. However, at the lower end of the photon spectrum the tail of the theoretical distri- 
bution is significantly narrower than the experimental one. Furthermore, at high photon 
energies, major discrepancies also appear because the contributions there correspond to 
processes where the outgoing neutron is either bound in the residual nucleus or it is in 
the continuum but carrying low kinetic energies. The first contributions have not been 
considered in our approach yet while the latter ones depend strongly on the details of the 
nuclear structure model. 

One may think that the situation could be improved by allowing for reasonable mod- 
ifications of the mean field potential. In Fig.|^ we show the results for dR^^^ /dq in '^''Ca 
obtained with the four different potentials defined in Table ^. They differ each other 
only in the proton or neutron depth of the Wood-Saxon well (parameter Vq in Eq. (|])). 
We see in the figure that none of the new potentials (B,C,D) give a photon distribution 
significantly wider than the A potential, used throughout this work. However, the be- 
havior above 110 MeV (low nuclear excitation energies) turns out to be very sensitive to 
the potential. The area below the curves depends also appreciably on the model and it 
increases when the neutron Wood-Saxon depth decreases. 

Now we focus our attention in the effects of the 2p2h correlations. Following the 
method described in sect. p73| , we include the effects of the ph self-energy in the medium, 
into the pure one-body SM differential rates of Fig. |l]. The results are shown in Fig. ^. 
We find an excellent agreement with the experimental distribution, specially in the low 
energy tail. Clearly, the inclusion of two body mechanisms has made the distribution 
wider and turns out to be essential to explain the tail of the experimental data. For 
comparison we also show the results for '^''Ca with the model of Ref. 0. In this reference. 
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q [MeV] 

Figure 2: Results for dR'-^ydq in ^^C, i^Q and ^OQa. Only one-body contributions have 
been included. Experimental data, taken from |0] and are given in arbitrary units. 
Our results have been normalized to the data at g = 105 MeV. In oxigen and calcium we 
have convoluted our differential decay rate with a Gaussian weight function of width 1.5 
MeV. 
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Figure 3: Results for dR^'^^^dq in '^°Ca with the different SM potentials defined in Table ^. 
The solid line corresponds to the potential A outlined in sect. ^]2| and used in the rest of 
this paper. Dash-dotted, dotted and dashed curves correspond to the potentials B,C and 
D respectively. 



the calculations were performed in nuclear matter and the LDA was used to obtain results 
in finite nuclei. There, though it was found a fair agreement for the integrated ratio R^'^'^ 
throughout the periodic table, it was pointed out the inability of the model to describe 
the tail of the photon distribution. In Ref. this problem was associated to the lack 
of high-momentum components in the FG nucleon wave functions. Here, we show that, 
in addition to the use of proper nucleon wave functions, one needs to incorporate 2p2h 
correlations to understand the experimental distribution. 

To discuss in more detail the effect of the inclusion of ph self-energy in the calculation 
of the differential decay rate, we show in Fig. ^ the results for dR^"'ydq in ^^O with and 
without the inclusion of 2p2h correlations. As it can be seen, these effects reduce the 
strength of the differential rate in the region of energies around the peak due to the pure 
one-body absorption, as it would do a FSI treatment with an optical potential. But, in 
addition these 2p2h correlations incorporate the contribution to the RFC process due to 
two-nucleons mechanisms in the low energy part of the photon spectrum. The net effect 
is a redistribution of the strength in the spectrum, with a variation of less than 10% in 
the integrated rate R^"'\ Similar effects have been also found in inclusive (e, e') processes 
at the quasielastic |]ll|, |jl9|, and A peaks ||2^ and also in inclusive {p,p') and {n,p) 



processes |TB[ . 

Besides, we would like to stress that the use of the same mean field potential for parti- 
cle and hole states in our calculations, guarantees the orthogonality of the ph states. It is 
common in the literature (|]T|, |TU[) to use different potentials for particles and holes (op- 
tical models for the particles and mean field potentials for the holes, etc..) and therefore 
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600 




q [MeV] 

Figure 4: Results for dR^^'^dq in ^^C, ^^O and '^^Ca. Two-body contributions have been 
now included. Experimental data, taken from and [Q, are given in arbitrary units. 
Our results have been normalized to the data at g = 105 MeV. For comparison in ^°Ca we 
have also shown the FG results presented in Fig. (3) of Ref. (thick line). The absolute 
value of the FG distribution around the peak is about a factor two or three greater than 
the SM one (thin line) and it has been normalized to the peak of the data. Additionally, 
the FG curve has been shifted 10 MeV to the left. 
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Figure 5: Results for dR^^^/dq in ^^O with (solid line) and without (dashed line) the 
inclusion of 2p2h correlation effects. The area below the curves are 9.2 x 10^^ and 9.8 x 10^^ 
respectively. 



the orthogonality between particle and holes states is lost. In these circumstances, there 
are new contributions to the transition matrix elements coming from the overlap between 
non-orthogonal wave functions. These new contributions, though not significant for high 
momentum transfer processes E^, can be important for processes like this where both 



energy and momentum transfered to the nucleus are small. 

The last aspect we want to discuss concerns the integrated ratio R^'^^ In Table ^ we 
present our results for this ratio and compare them with experiment and with those of 
the FG model of Ref. §. 

The results obtained within the framework described up to here (column labeled "Con- 
tinuum" in the table) turn out to be around a factor of two smaller than the experimental 
ones and those of Ref. 0. However, it is worth to point out here that, as shown in Table ^, 
the integrated rate is rather sensitive to the details of the SM nuclear potential used, and 
that a small variation of the model could account for such discrepancies. Apart from this 
point, other possible reasons for the low rates found are the following: 



1. We have only considered the Kroll-Ruderman term for the elementary process 
7T~p 7n. The momentum of the pion on its bound orbit though small is not 
strictly zero, and one should add to the dominant Kroll-Ruderman term new am- 
plitudes proportional to the pion momentum which could increase the integrate 
RFC rate. In any case these new contributions are not expected to be important 
enough to explain the disagreement. 
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[%] 




Ref. 11 


Continuum 


Discrete Total 


Experiment 


12q 


1.19±0.16 


0.62 


1.17 1.79 


1.84 ± 0.08, 1.92 ± 0.19, 1.6 ± 0.1 


160 


1.53 ±0.08 


0.92 


0.80 1.72 


2.27 ±0.24, 2.24 ±0.48 


4"Ca 


1.91 ±0.29 


1.17 


0.87 2.04 


1.82 ±0.05 



Table 4: Results (in percent) for R^'^\q) in ^^C, ^^O and ^°Ca. The second column contains 
the results obtained with the FG model of Ref. Our results are given in the next three 
columns: third column corresponds to the area below the curves in Fig. ^ fourth column 
gives the contribution of the RFC processes in which the final neutron is bound in the 
residual nucleus and finally the fifth column is the total rate, sum of the two previous 
values. Experimental values are taken from references: [^, [0 and for ^^C, and 
for 1^0, and M for ^^Ca. 



2. The determination of the weights Uni (see Table 0) may be subject to theoretical 
uncertainties. It would be interesting to disentangle experimentally the capture 
from different pionic atomic states to allow a direct comparison with the theory, 
free of the assumptions made on the values of Uni (for more details see Ref. 0). 

3. Pion capture processes where the outgoing neutron is bound in the residual nucleus 
instead of going to the continuum lead to Dirac's delta peaks in the photon spectrum 
which also contribute to the integrated ratio R^"''^ and have not been considered yet. 
Indeed, one can easily evaluate these contributions by taking into account that the 
integral over in Eq. (|12D is now a discrete sum over the unoccupied single neutron 
states (see Ref. |15|) and the radial wave function Rp in Eq. (^) corresponds now 



to the solution of the radial Schrodinger equation for a negative energy ep. 

In the fourth column of Table ^ we give the contributions from these discrete states 
to the integrated rate. We see that they are of the same order of magnitude than 
the continuum ones and lead to a reasonable agreement with the experimental data. 
To describe properly the magnitude and the position of each individual discrete 
transition, one would need a more sophisticated model for the residual ^^B, ^^N 
and nuclei than the one used here, as it was the case when we discussed the 
continuum integrated rates. However, the sum of all the discrete and continuum 
contributions is rather model independent as it can be seen in Table |^. Results 
with potentials B, C and D agree within around 4%, while the value obtained 
with potential A is well below the previous ones. This is because the proton well 
for potential A is more attractive than for the rest of potentials, and therefore the 
overlap between the holes and pion wave functions is smaller. Changes in the proton 
potential modifies both the continuum and discrete contributions, and then the total 
rate. However, variations of the neutron potential depth, leaving fixed the proton 
one, make some single-particle states to be bound or not. In this last case, low 
energy resonances appear in the spectrum. Thus, the net effect is an exchange of 
strength between the discrete and continuum contributions to the RFC rate, leaving 
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Shell-model 



Potential /Density 


Continuum 


Discrete 


Total 


LDA Fermi gas 


A 


1.17 (1.25) 


0.87 (0.92) 


2.04 (2.17) 


2.20 


B 


1.44 (1.54) 


0.97 (1.03) 


2.41 (2.57) 


2.62 


C 


1.76 (1.87) 


0.57 (0.61) 


2.33 (2.48) 


2.59 


D 


1.14 (1.21) 


1.28 (1.36) 


2.42 (2.57) 


2.66 


2pF 








2.74 



Table 5: Comparison of the results (in percent) for R^"'\q) in '*°Ca for the same potentials 
of the Table ^ and the FG model of Ref. ||^. In the first column we give either the 
Wood-Saxon potential used to compute the SM results or the proton and neutron center 
densities needed as input in the FG model. Densities A-D refer to those obtained from the 
corresponding Wood-Saxon potentials A-D and 2pF to the two parameter Fermi density 



used in the Ref. P| and taken from Ref. ||3^. The 2pF value does not match the one 



shown in Table |^ because here RPA correlations have not been included and the error 
bars due to uncertainties in the strong absorption widths, F";"*, have been dropped out. 
To compare with the FG model we also show in brackets the results of the SM calculation 
without taking into account 2p2h correlations. 



almost unchanged the total value. 

From the above discussion we see that reasonable changes in the proton potential 
lead to about 15-20% uncertainty in the SM results. Moreover, one gets an addi- 
tional 4% uncertainty associated to the neutron potential, thus we conclude that the 
SM gives integrated RFC rates with an uncertainty of about 20-25%. Similar results 
are also found for carbon and oxigen. This uncertainty is basically due to two facts: 
i) the RFC process takes place at the nuclear surface and ii) the main contribu- 
tion to the integrated RFC rate comes from nuclear excitation energies rather low 
(?a 5 — 10 MeV). For processes sensitive to the whole nuclear volume and involving 
larger excitation energies (see for instance calculations at the quasielastic peak for 
inclusive (e, e') scattering |]1T|) the SM results are much more independent of the 
details of the mean field. 

4. We have not included RPA correlations in our model. In Ref. RPA lead to a 
reduction of about 40% in the integrated ratios without changing the shape of the 
photon distribution. In Table ^ we compare the SM results, without including 2p2h 
correlations (in brackets), to those obtained with the LDA FG model of Ref. 
without including RPA correlations (last row). We see that the FG model value 
(2.74) is significantly higher than the SM ones, part of this discrepancy being due 
to the use of different proton and neutron center distributions in both types of 
calculations. To make more reliable the comparison, we also show the LDA FG 
results obtained using the proton and neutron densities derived from the different 
Wood-Saxon potentials used in the SM calculations. We can see now that the FG 
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model gives similar results to those of the SM and it is also subject to the same type 
of uncertainties (20-25%) associated to the precise details of the nuclear dynamics. 
Similar results are obtained in carbon and oxygen. 

However, we would like to point out that the LDA FG results are in most cases above 
the SM ones. This can be due to the undoubted deficiencies of the FG picture of 
the nucleus to describe processes where low nuclear excitation energies are involved. 
But one could also try to explain this disagreement (or part of it) by questioning 
the validity of the LDA to describe a surface process like the one studied here. In 
any case, because of the zero range character of the Kroll-Ruderman interaction 
considered here, LDA should work much better than what one could expect. 

In any case, and in view of the results in Table |], one might think that the effect 
of the inclusion of the RPA correlations could be similar in both models. The fact 
that the SM result shows an uncertainty comparable to the expected RPA effects 
prevents us to try to include such correlations before having a more reliable nuclear 
model. 



4 Summary and conclusions 

In this paper we have reanalyzed the RFC process in light and medium nuclei. We have 
presented a simple continuum SM which describes successfully the shape of the photon 
distribution. We have paid a special attention to the study of the low energy part of the 
photon spectrum which could not be fully understood within the different models found 
in the literature and in special in the model of Ref. 0, despite of being the only one 
able to describe the experimental integrated ratios R^"'^ thorough the periodic table. We 
have found that the effect of 2p2h correlations play an essential role to describe the tail 
at low energies of the photon distribution and that these effects together with the use of 
nucleon wave functions with non-vanishing high momentum components allow us for a 
correct description of this part of the spectrum. 

We have also found that the region of high photon energies depends strongly on the 
details of the nuclear model used. In that region when the outgoing neutron is bound, the 
photon spectrum exhibits peaks characteristic of the energy spectrum of the final nucleus. 
Also in this region, when the emitted neutron goes to the continuum, there are peaks 
corresponding to resonances of the final residual nucleus. Much effort have been invested 
in the past to extract information about these nuclear bound and resonances states pi, pi 



and ||I0|. We think that a first necessary step to attack this interesting problem is a correct 
understanding of the quasifree background, studied here, which requires the inclusion of 
two body mechanisms and to deal with orthogonal particle-hole states. Thus, results 
presented in this paper could be used in future to isolated the contribution of such nuclear 
states from the experimental data. This would clearly represent an improvement with 
respect previous works where the subtraction of the quasifree background was performed 
by means of optical model calculations or the non-theoretical founded pole model (see 
Ref. m). 
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Finally, we have obtained values for the integrated ratio R^^^ in good agreement with 
those found with the FG model of Ref . . We have also discussed different mechanisms 
which could account for the remaining discrepancies with the experimental data. None 
of them are expected to make narrower the photon distribution because, as we have 
mentioned above, the tail of the distribution is mainly due to 2p2h effects and high 
momentum components in the nucleon wave functions and these physical effects will be 
present in all new contributions to the RFC rate. 
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Appendix 

Matrix elements of the electromagnetic operators 

In this appendix we give the matrix elements needed to calculate the width of the RFC 
process as given by Eq. (H). As it can be seen in Eq. (||), J(r) is in our case a one-body 
operator in the nuclear space. Thus both Tj^ and Tj^^ operators are also. We define 
the auxiliary operators U^^^"^^^ by integrating the operators T*^''™^^ over the pionic degrees 
of freedom and projecting the result onto a well defined angular momentum basis. 



iel,mag/ 



J'M' 



(JMJ.M.I J'M%(0|rj'^^^(g)|7r-; J,M,) (19) 

MMi 



where the pion states [vr; Jj) and |0),r where defined in section ^]2[ The Wjj™^|/j,(g) oper- 
ators act now only on nuclear states and have rank J'M'. The reduced matrix elements 
of the T*^''™*^^ operators can be related to those of the W^^'^^s operators, 

{ph-\jf\\Tf'-^{q)U) = {-iy^-'-'HpmtjJ{q)\\h), (20) 

where the particle and hole states with third isospin component tp and th respectively are 
given by: 

1 1 

\P) = Rpir,ep)\-lp;jpmp)\-,tp) (21) 
\h) = Rh{r,eh)\hh;jhmh)\^,th). (22) 
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As we discussed in sect. pl2| , the radial wave functions Rp and Rh are obtained by solving 
the Schrodinger equation with the potential described in that section, for positive (ep) 
and negative (et) energies, respectively. The values of eh are determined by the quantum 
numbers Uh, h and jh of the occupied shells in the parent nucleus, whereas the positive 
energies are determined by energy conservation. The particle radial wave functions Rp 
are normalized to the Dirac's delta of energies. 

By using Eqs. (0)-(||) and Eq. ( p!OD one can find the following expressions for the 
single-particle Wfj^jFj^, (q) operators: 



- V2£(l + |l|)^*_ (24) 
1 



x-j E J + ^s,-iJj+siqr)RljSr) [[Yj+sir) (g) a]j 0YjXr)] 

•J s=±l 



JfMf 



where r is now the nucleon spatial coordinate and J = y/2J + 1. 

By using Racah algebra to compute the reduced matrix elements of the operators in 
Eqs. (H-H) we get finally: 



X3p3jjjfi{ip + h + j + j.;)[^^l \_ 

Kp + {-l)lv+lH+Jf+l^^^ I J J J\ 

^ ^ 1-10 



JAJf + 1) 

/ drr^Rp{r)3.j{qr)R^{r)Rh{r) (25) 

JO 



Kp + {-iyr+iH+jf+i^^ / J \ 

1-10 



j^^jf{jj + i) 

POO 

X / drr^Rp{r)[J]j+i{qr) - {J + l)jj^i{qr)]R^{r)Rh{r) 
Jo 
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where we have defined k, 
even (odd) values of n. 
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